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In  this  work  we  Investigate  numerically  the  problem  of  the  source 


excitation  of  an  open,  parallel-plate  waveguide.  The  following  assump- 
tions are  made  for  the  source  current  1)  the  current  Is  oriented  In  the 
y-dlrectlon,  2)  It  Is  located  at  x » 0,  3)  j|ihere  Is  no  variation  In  the 
y-dlrectlon.  A)  and  the  current  has  exp  (kSz)  behavior  along  the  longi- 
tudinal z-dlrectlon.  We  provide  graphical  output  for  the  EM-fleld 
components  as  functions  of  a longitudinal  propagation  constant  and 
transverse  coordinates  and  then  discuss  these  results. 
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I.  INTRODUCTION 


In  the  previous  report  [1],  we  derived  analytical  expressions  for 
the  source  excitation  of  an  open  parallel-plate  waveguide.  However, 
these  formulas  were  very  complicated,  and  it  became  necessary  to  evaluate 
them  numerically.  The  purpose  of  this  report  is  to  present  the  numerical 
results.  The  computer  program  contained  in  Appendix  A was  written  and 
used  to  obtain  the  field  distribution  as  a function  of  the  longitudinal 
propagation  constant  and  the  transverse  coordinates.  The  numerical  out- 
puts are  presented  in  graphical  forms.  The  Cyber  175  at  the  University 
of  Illinois  was  used  for  all  of  the  numerical  studies. 

The  organization  of  the  report  is  as  follows:  Section  II. contains 
a statement  of  the  problem  and  the  basic  formulation.  Section  III 
presents  the  real  and  Imaginary  parts  and  the  amplitude  of  the  component 
field  distribution  as  functions  of  several  parameters  In  graphical  form 
and  a detailed  discussion  of  the  numerical  results.  Finally,  Section  IV 


Is  the  conclusion. 


II.  STATEMENT  OF  THE  PROBLEM  AND  BASIC  FORMULATION 

In  this  section  the  fields  due  to  a vertical  current  located  inside 
an  open,  finite  waveguide  are  investigated.  The  geometry  of  the  problem 
considered  is  shown  in  Figure  1.  This  structure  consists  of  two  perfectly 
conducting  plates  with  separation  2H  located  in  a homogeneous  and  Isotropic 
medium.  A Cartesian  coordinate  system  with  its  y-axls  normal  to  the  plates 
is  erected.  Both  plates  are  infinite  in  the  z-dlrection  and  finite  in  the 
x-dlrectlon  with  length  2L  as  shown  in  Figure  1.  All  figures  appear  at  the 
end  of  Chapter  II.  The  current  is  oriented  in  the  y-dlrectlon  and  is 
defined  as 

J ■ y6(x)  exp  (IBz)  , (1) 

where  6 is  the  propagation  constant  in  the  z~direction,  and  6(x)  is  the 
delta  function.  In  [1]  using  the  vector-potential  approach  and  the  Wlener- 
Hopf  technique,  we  obtained  a solution  for  the  problem  at  hand  in  a general 
form  for  any  parameters  with  one  restriction:  kL  must  be  much  greater 
than  1,  l.e. , 

kL  » 1 , (2) 

where 


and  e,  u are  the  homogeneous  media  parameters.  Using  the  solution  which 
was  obtained  in  [1],  we  will  perform  a numerical  calculation  for  the  case: 

W - ^ - 0.16670  , (4) 

^ - 5 , (5) 

where 

Xo  - (6) 


jv 


2 


Is  the  free'-space  wave  length.  Because  of  the  limitations  of  Equations 
(2),  (3),  and  (5)  we  calculated  numerical  results  for 


0 < e < 0.93  , (7) 

**  6 
where  6 Is  the  normalized  propagation  constant  3 * . It  can  be 

readily  proved  from  Equations  (3) , (4) , and  (5)  that  If  3 la  In  the  region 
[3q,  0.93] , where 

3g  - 0.80008997  , (8) 


then  only  the  first  mode  can  propagate  In  the  x-dlrection;  therefore,  the 

electromagnetic  field  has  only  three  components  [E  , H , H ] . When  3 

y X z 

Intersects  the  point  3q  and  goes  to  zero,  the  second  mode  begins  to 
propagate,  and  the  EH-fleld  consists  of  five  components  [E  , E , E , H , H ] . 

X Z X z 

Rewriting  the  field  solution  from  [1]  for  the  case,  when  no  more  chan  two 
modes  can  propagate,  we  arrive  at  the  EM  field: 


E'(x,y,z)  ■ E(x,y)  • exp  (13z)  (9) 

H'(x,y,z)  - H(x,y)  • exp  (i3z)  (10) 

A ^ A 

E(x,y)  - xE^  + yE^  + zE,  (11) 

X y z 

H(x,y)  - xH^  + (12) 


4 


it 


I 

1 


i 

i 


• - p ] ^1  *°®  ( ^ ^ ^2  ^ * h)  * * l) 


y'ep 


^exp  |la 


\'N 


a 


■ lyaFj^  sin  |a^j  + iSFj  cos  |Tr^]  sin  (s  * ^|  - 0*5  * exp  [ ia  ^ • j-^j 


(1 


[M^^(k)]^  exp  (i2a)  r 2b[M^_j,(oi^)  exp  . (1 

‘ )a  r ' 

2bM^^(k)M^^(o^)  exp 


1 (1  + T^)a 


ia  11  + ^ 


a • Q 


(1 


2 — /.»%a^.  V iw  L V 2“\ 


(1  + • 


|[M^^(a^)]"  exp  (12  fa^) 


2o, 


ATj^b  (Oj^  + b)' 


.1  + Tj^  “l  J 


1 . (2 


k)  - (a,+b) 


j^+b)  • exp  ^ 


b(2-C  + In(^)  + i^)  - ^ + I (^  - arc  sin  -) 

b 2 2 ni;2  " 


(2 


••  i/2  • OL^  • 


1 - ia. 
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* I 


^ exp^  i|a^(2-c  + In  (I")  ■»•  iy) 


n-2 


xrrijj)  ’ 


(2 


(2: 


a ■ kL  ■ lOir 


S 12 


(2 


exp  (i2a)  j^l  + ^ 


exp  (-  i -p 


It  should  be  mentioned  that  we  Investigated  the  lossless  medium  case; 
therefore.  In  the  region  Sq  < 8 < 0.93,  the  propagation  constant  for 
the  second  mode  has  only  an  Imaginary  part.  Because  we  neglect  terms 
which  decrease  exponentially,  our  results  for  the  above  mentioned  region 
of  8 reduce  to: 

1)  E - E - 0 and 
X z 

2}  more  simple  expressions  for  the  other  three  components  of  the  field. 
We  apply  numerical  analysis  only  over  the  regions  0£y<H,0<x<L. 

For  the  remainder  of  the  waveguide,  one  can  obtain  results  using  the 
correlations : 


- E^(-x,y) 

; E^(x,y)  - 

- E^(x,-y) 

Ey(-x,y) 

; Ey(x,y)  - 

Ey(x,-y) 

; E^(x,y)  - 

- E^(x,-y) 

H^(-x,y) 

; H^(x,y)  - 

H^(x,-y) 

- H^(-x,y) 

; H^(x,y)  - 

It  Is  Interesting  to  note  that  E , H are  continuous  and  that  E , E , H 

jC  z y z X 

are  discontinuous  when  8 crosses  8^  (or  more  exactly:  they  are  exponentially 
decreasing).  It  Is  also  of  Interest  to  determine  the  character  of  the 
behavior  of  the  x-component  of  Poyntlng's  vector.  As  one  can  easily  see 
from  the  previous  expressions  for  the  EM  fields,  the  x-component  of 
Poyntlng's  vector  for  the  second  mode  Is  proportional  to  and  goes  to 
zero  when  8 -*•  8g.  From  this,  one  finds  that  the  energy  flow  In  the 
x-dlrectlon  Is  continuous  when  8 Intersects  the  point 
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Figure  2.  The  real  and  Imaginary  parts  of  an  x-component  of  the  electric  field  as  functions 
of  a longitudinal  propagation  constant  for  points  of  view:  -7  = 0.1;  — - 0.0,  0 
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Figure  5.  The  real  and  imaginary  parts  of  an  x-component  of  the  electric  field  as  functions 
of  a longitudinal  propagation  constant  for  points  of  view:  ^ = 0.9;  ^ = 0.0,  0 
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eal  and  imaginary  parts  of  a y-component  of  the  electric 
longitudinal  propagation  constant  for  points  of  view:  — 
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eal  and  imaginary  parts  of  a y-component  of  the  electric  field  as„functlons 
longitudinal  propagation  constant  for  points  of  view:  ^ = 0.6;  = 0.0,  0 
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Figure  11.  The  real  and  imaginary  parts  of  a z-component  of  the  electric  field  asYfunctlons 
of  a longitudinal  propagation  constant  for  points  of  view:  ^ = 0.4;  — = 0.0,  0 
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Figure  19.  The  real  and  imaginary  pares  of  a z-component  of  the  augneCic  field  as^functions 
of  a longitudinal  propagation  constant  for  points  of  view:  -r  • 0.4;  „ ■ 0.0,  0 
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eal  and  imaginary  parts  of  a z-componenC  of  the  magnetic  field  as^functlons 
longitudinal  propagation  constant  for  points  of  view:  — 0.9;  ^ » 0.0,  0 


le  moduleof  a y-oomponent  of  the  electric  field  as  a function  of 
ngitiidinal  propagation  constant  for  points  of  view:  y “ 0.6,  0 


The  nodule  of  a z-conponent  of  the  magnetic  field  as  a function  of 
longitudinal  propagation  constant  for  points  of  view:  ^ = 0.6,  0 


gure  29.  The  real  and  imaginary  parts  of  a^y-component  of  the  electric  field  as  functions 
of  an  x-coordinate  for  Z = 0.75;  — = 0.0,  0.5. 


Figure  32.  The  real  and  Imaginary  parts  of  a^z-component  of  the  magnetic  field  as  functions 
of  an  x-coordlnate  for  ' = 0.4;  — = 0.0,  0.5, 
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III.  NUMERICAL  STUDY  OF  THE  PROBLEM 


Because  the  analytical  expressions  derived  In  the  second  section 
are  rather  complicated  and  difficult  to  analyze,  we  numerically  evaluated 
the  solutions  using  a digital  computer.  The  results  are  presented  In  this 
section  In  graphical  form.  This  section  consists  of  two  parts:  In  the 
first  we  discuss  the  field  components  as  functions  of  the  longitudinal 
propagation  constant;  In  the  second  > as  functions  of  the  transverse 
coordinates. 

A.  Real  and  Imaginary  Parts  and  Amplitudes  of  the  Field  Components 


As  Functions  of  the  Longitudinal  Propagation  Constant 

The  graphics  that  are  supplied  In  this  section  were  plotted  with  a 
step  for  B equal  to  0.005.  The  point,  B ■ 0.80009,  also  was  used.  The 
figures  were  plotted  using  188  points.  The  output  for  the  real  and 


Imaginary  parts  and  the  amplitudes  of  all  five  field  components  for  eight 


X Y 

observation  points:  — ■ 0.1,  0.4,  0.6,  0.9  and  — ■ 0.0,  0.5  are  presented. 


From  the  figures  It  Is  observed  that  E^,  H^,  are  dominant  components. 


In  Figures  6-9  the  real  and  Imaginary  parts  of  are  given.  By  comparing 


Y Y 

these  results,  one  can  observe  for  points  — - 0.0  and  — - 0.5  that  the 


E , H , H - field  behaviors  as  functions  of  B are  more  complicated  In 
y X z 


the  middle  of  the  waveguide.  This  Is  hardly  surprising  In  view  of  the 
largest  contribution  of  the  second  mode  for  the  above-mentioned  com- 


ponents for  ^ ■ 0.0.  When  the  point  of  observation  approaches  - 0.5, 


the  contribution  of  the  second  mode  for  those  components  decreases  and 
has  a limiting  value  equal  to  zero.  In  Figures  6-9,  14-17,  22-25  we 


observe  the  step  changing  for  the  E^,  components  at  B ■ Bq.  The 
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second  mode  Is  responsible  for  this  misbehavior.  When  ^ 0.5,  Che 

n 

concribuClon  of  Chat  step  decreases,  and  for  ^ - 0.5,  it  equals  zero. 

One  can  observe  that  Che  curves  are  smooch  ac  S * Sq.  Ic  should  be 
mencloned  chac: 

(a)  The  dlsconClnuous  behavior  Is  observed  because  we  have 
neglecCed  Che  exponentially  decreasing  cerms; 

(b)  In  a lossy  medium,  Che  rate  of  decrease  would  be  less.  The 

H^-componenC  (see  Figures  18-21,  26)  doesn't  have  the  step  behavior 
because  the  contribution  of  Che  second  mode  Is  proportional  Co  a^, 
which  goes  to  0,  when  8 8^.  The  other  two  components,  and  E^,  are 

equal  to  zero  on  the  x-axls  and  have  their  largest  contribution  for 

^ ■ 0.5  (see  Figures  2-5,  10-13).  E Is  smooth  and  E exhibits  Che  step 

H X Z 

behavior  at  8q.  The  figures  demonstrate  that  the  complexity  of  Che 
curves  occurs  approximately  In  the  region  0.7  < 8 < 0.8,  where  we 
observe  a sharp  peak,  which  Is  due  to  the  resonance.  In  the  region 
0 < 8 < 0.5  Che  figures  demonstrate  Che  very  smooch  character  of  Che 
curves . 

B.  Real  and  Imaginary  Parts  of  the  Field  Components  as  Functions 
of  Transverse  Coordiaatea 

In  this  section  we  present  the  real  and  Imaginary  parts  of  the 
dominant  field  components  E^,  as  functions  of  the  x-coordlnate 

for  two  values  of  ^ ■ 0.0;  0.5  and  two  values  of  8 ■ 0.4;  0.75.  The 
graphic  output,  shown  In  Figures  28-33,  was  obtained  using  Che  results 
of  calculations  for  81  points  of  ^ In  the  region  [0.1,  0.9]  (step  ••  0.010). 
All  the  graphs  have  very  smooch  characteristics.  For  8 ■ 0.75  they  have 
slightly  mors  complicated  form  Chan  for  8 *0.4.  As  mencloned  In  Che 
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I 

I 

i 

previous  section,  for  ^ ■ 0.5  the  contribution  of  the  second  mode  equals  j 

zero.  We  see  that  the  amplitudes  of  the  curves  are  constant  for  the  j 

entire  region  of  view.  For  the  ^ - 0.0,  the  field  components  are  sums  of 
the  contributions  of  two  modes.  One  can  observe  that  the  amplitudes 
of  the  cur'/es  are  changing  along  the  x-dlrectlon. 


IV.  CONCLUSIONS 


In  chls  report  the  problem  of  a source  excitation  of  an  open  I 

parallel-plate  waveguide  was  developed.  Extensive  numerical  results  for 
the  field  components  In  the  waveguide  as  functions  of  several  parameters 
of  the  waveguide  and  propagation  constant  were  supplied. 
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APPENDIX 

SOURCE  EXCITATION  OF  AN  OPEN,  PARALLEL-PLATE  WAVEGUIDE  PROGRAM 

A complete  program  for  source  excitation  of  an  open,  parallel-plate  wave- 
guide program  Is  presented.  The  computer  program  provides  three- 
dimensional  data-storage  for  the  real  and  Imaginary  parts  of  five  components 
of  the  EM  field.  Data  were  obtained  for  ^ between  0.1  - 0.9  with  step  0.1; 

between  0.0  - 0.9  with  step  0.1;  and  S - propagation  constant  between 
0.0  - 0.93  with  step  0.005  plus  (•)  0.80009.  These  data  were  used  to 
plot  EM  field  components  as  functions  of  the  propagation  constant.  The 
program  can  be  readily  modified  to  obtain  data  for  plotting  the  EM-fleld 
component  as  a function  of  the  x-coordlnate . 
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PROGRAM  AFIELD(INPUT»QUTPUT»TAPE3»TAPE1=INPUT) 

COMPLEX  BBK » BBAL » T1 1 FI » F2  »EXO»EYO»EZO»HXO»HZO»  CONST  »AA»EiB 
DIMENSION  XX(9) f YY(IO) »BETAB< 188) » REXO ( 10 r 9 » 138 ) » AMEXO < 10 r 9 , 188 > » 
*REY0(10»9»188) »AMEY0(10»9»188) ,REZ0<10»9»188) »AMEZ0(10»9»188) » 
*RHX0<10r9»188) »AMHX0(10f9»188) rRHZO ( 1 0 r 9 r 188 ) »AMHZO( 10r9f 188) 

RE AD  < 1 » 2 ) W » YBR I » YBRS  f YBRF  f XBRI » XBRS » XBRF 
2 F0RMAT<F7.5»6(F7.3)) 

M=1000 

PI=3. 141592654 
C0NST=CMPLX<0. >1.  ) 

C0N=2. 11593152 
EETAB<1)=0. 

DO  70  I=l»160 
BETAB ( 1+ 1 ) =BETAB ( I ) +0 . 005 
70  CONTINUE 

BETAB(162)=0. 80009 
BETAB< 163)=0.805 
DO  75  1=1 »25 

75  BETAB(I+163)=BETAB(I+162)+0.005 
CONTINUE 
DO  98  N=l»161 
BETA=BETAB(N) 

K=N 

A=10.*PI*SQRT<1.-BETA**2) 

B=A#U/PI 

ALFAl=SaRT<B**2-l. ) 

DD=PI*ALFA1/W 

CALL  BEl(AfB>BBK>PI>MfCONf ALFAl) 

CALL  BE2(Af BrBBAL>PIfMfCON>ALFAl ) 

CALL  TE1(A»B»T1»BBK»PI) 

CALL  F 12 ( A r B » T1 r ALFAl » BBK » BBAL » F 1 » F2 ) 

X=XBRI 

J=1 

30  Y=YBRI 
1 = 1 

20  EXO=(PI*ALFAl/<  10.)KW**2)  )*F2)(cSIN(PI*Y)*SIN(DD)tcX) 

EY0=10.*PI*<F1*C0S<A*X)  + <1  ♦-(  1 ,/( 100.»U#*2)  ) ) !lcF2*C0S <PI*Y ) * 
♦C0S(DD#X)-CEXP<CMPLX<0. »A»X) )/(2.*A) ) 
EZ0=-C0NST*BETA*<PI/U)#F2*SIN<PI*Y)*C0S(DD*X) 

HX0=-BETA*10.*PI*(Fl*C0S(A»X)+F2#C0S(PI»Y)*C0S(DD*X)-( 1 ./<2.*A> )» 
*CEXP(CMPLX(0. »A*X) ) ) 

HZ0=C0NST»A*F1*SIN< A*X ) +C0NST*DD*F2*C0S (PI*Y ) *SIN  < DD*X ) - 
♦ 0 . 5*CEXP ( CMPLX ( 0 . » A*X ) ) 

REXO ( I » J f K ) =REAL  < EXO ) 

AMEXO(I»J,K)*AIMAG<EXO) 

REYO(I» JrK)=REAL(EYO) 


f 
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AMEYO< I» J.K)=AIMAG(EYO) 
REZO<I»JfK)=REAL(EZO) 
AMEZO< I» JfK)=AIhAG<EZO) 
RHXO  < I . J . K ) =REAL ( HXO ) 
AMHXO< If J»K)=AIMAG<HXO) 
RHZO  < I f J f K ) =REAL ( HZO ) 
AMHZO< If JfK)=AIMAG<HZO) 
Y=Y+YBRS 
1 = 1 + 1 


IF(Y.LE.YBRF) 

GO 

TO 

20 

X=X+XBRS 

J=J+1 

IF(X.LE.XBRF) 

GO 

TO 

30 

98  CONTINUE  1 

DO  99  L=162fl88  j 

BETA=BETAB<L) 

K=L 

A=10.*PI)KSaRT(l  .-BETA*#2) 

X=XBRI 

J=1 

B=A*U/PI 

CALL  BE(AfBf BBfPIfM)  ; 

X=XBRI 

17  AA=CEXP(CMPLX(0. f A#X) ) 

EY0=-0 . 5* < AA-BB*COS  ( A»X ) ) /SORT  ( 1 , -BETA)tt*2 ) 

HXO=-BETA*EYO 

HZ0=0.5*<-SIGN(1. fX)»AA+CONST*BB#SIN<A*X)  ) 

DO  18  I=lflO  ' 

REXOdf  JfK)=0. 

AMEXO<If JfK)=0. 

REZOdf  JfK)=0. 

AMEZOdf  JfK)=0. 

REYOdf  JfK)=REAL(EYO)  j 

AMEYOdf  JfK)=AIMAG(EYO) 

RHXO  < I f J f K ) =REAL ( HXO ) 

AMHXOdf  JfK)=AIMAG<HXO)  ' 

RHZO (I f J f K ) =REAL ( HZO ) 

18  AMHZOdf  JfK)=AIMAG(HZO)  i 

X=X+XBRS 

J=J+1  j 

IF<X.LE.XBRF)  GO  TO  17  i 

99  CONTINUE  : 

URITE(3f 101)  REXOfAMEXOfREYOf AMEYOf REZO fAMEZOf RHXO fAMHXOf RHZO f AM  1 

♦ HZO  I 

101  FORMATdOFS.S)  1 

DY=0.1 

DX=0,1  I 

YY(1)=0. 

DO  50  I=lf9 

YYd  + l)=YY<I)+DY  I 

50  CONTINUE  ’ 
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XX(1)=0.1 
DO  60  1=1 »8 
XX(I+1)=XX<I)+DX 
60  CONTINUE 

WRITE(3»101)  YYrXXrBETAB 

STOP 

END 

SUBROUTINE  BE 1 ( A f B » BBK » PI f H » CON f ALFA 1 ) 

COMPLEX  FfBBK 
AM1=0. 

DO  10  I=2fM 
AM1=AM1+B/I-ASIN(B/I) 

10  CONTINUE 

F=CEXP(CMPLX(-B*PI/2. fB*(CON-ALOG(B) )-PI/2.+A+AMl) ) 

BBK=(ALFA1+B)*F 

RETURN 

END 

SUBROUTINE  BE2< A»Bf BBALrPI »M»C0N»ALFA1) 

COMPLEX  FfBBALfD 
AM1=0. 

DO  10  I=2»M 

AM1=AM1+ALFA1/I-ASIN<ALFA1/SQRT( I**2-l . ) ) 

10  CONTINUE 

F=CEXP<CMPLX<-ALFAl!(tPI/2.  »ALFA1*(C0N-AL0G<B)  )+AMl+ALFAl*A/B)  ) 
D=CMPLX<1. f-ALFAl) 

BBAL=F*D)(tSQRT<2.  )*ALFA1/B 

RETURN 

END 

SUBROUTINE  TEl ( Af BrTl »BBK»PI ) 

COMPLEX  TlfBBK»Dl 
D=B*SQRT(PI)/SQRT(A*2. ) 

D1=CMPLX<1.+Df-D) 

Tl=BBKJtC!K2*Dl 

RETURN 

END 

SUBROUTINE  F12< Ar B » T1 . ALFAl f BBK f BBAL » FI , F2 ) 

COMPLEX  Tl»BBK»BBAL»Fl»F2fDl 

Dl=ALFAl)(t<l.+Tl)»<BBALJ«c»2*(4.*Tl*B/(l.+Tl)-<ALFAl+B)**2/ALFAl  )/ 
♦<2,*ALFA1)-1. ) 

Fl  = (BBK*#2/( < l.+Tl)*A) )*(1,+2.*B*BBAL**2/D1 ) 

F2=2 . *B)KBBK»BBAL/ < A#D1 ) 

RETURN 

END 

SUBROUTINE  BEIAr Br BBf PI f M) 

COMPLEX  TfF»BB 
C0NST=1 .115931516 
AM1=0. 

DO  10  I=1»M 
AM1=AM1+B/I-ASIN(B/I) 

10  CONTINUE 

T=CEXP<CMPLX(-BiltPIf2.*(B*<C0NST-AL0G(B)  >+A+AMl)  ) ) 
F1=SQRT<PI)*B/SQRT<2.*A) 

F=CMPLX(1.+F1»-F1) 

BB=2.*T/(1.+T*F) 

RETURN 

END 
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